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Patent Application 

FOR 

METHOD FOR CALCULATING AND DISPLAYING THE ISOTHERMAL 

CONTOURS PRODUCED BY A LASER 

Related Application 

This application is a divisional application of, and claims the benefit of, 
copending application U.S. Serial Number 10/271,706, filed on October 15, 2002, which 
claimed priority to U.S. Provisional Application No. 60/329,769, filed on October 16, 
2001, and U.S. Provisional Application No. 60/366,774, filed on March 22, 2002, both 
of which were incorporated in their entirety by reference. 

Field of the Invention 

The present invention relates to a method for calculating and displaying the 
isothermal contours produced by a laser in a sample. The method can be carried out 
using the described microscope and a computing apparatus for performing the 
calculations and displaying the isothermal contours. 

Background of the Invention 

Recent advances in biology and medicine have led to the development of laser 
beam microsurgery on cells. The laser beam is well adapted to micro-manipulation of 
small objects, such as single cells or organelles. It provides the advantage of non- 
contact ablation, volatilization, steriUzation and denaturing, cutting, and other forms of 
thermal and actinic-light treatment. The four parameters of focal spot size, laser 
wavelength, pulse duration, and laser power, provide a variety of regimes suitable for 

r 

different appHcations. One example use of laser beam microsurgery is the application of 
laser beams to the treatment of the mammalian oocyte and embryo. However, laser 
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beam microsurgery in a nurnber of inverted or upright microscopes can be utilized for 
many different surgical or medical applications. 

In further detail of the example application, the early-stage mammalian embryo 
5 is contained within a protective layer, the zona pellucida ("ZP"). , The ZP is relatively 
analogous to the shell of a hen's egg. This proteinaceous ZP layer is of varying 
thickness, typically 10 to 20 |jim, and of varying hardness. The embryo remains within 
the ZP during development from the single-cell to the blastocyst stage, at which point 
the embryo breaks out of the ZP and implants itself into the uterine wall. 

10 

It has been found that certain embryos, typically those from older niothers or 
embryos that have been frozen for storage, frequently have much tougher ZP layers than 
younger-origin or untreated embryos. Consequently, when the time comes for the 
embryo to emerge from the ZP, there may be a significant impediment in the tougher 
15 layers, which have to be traversed. If the embryo fails to hatch in the limited time 
available, it will be lost and fertility will fail. 

Assisted hatching derives from the observation that fertility can be augmented by 
generating holes or gaps in the.ZP through which the embryo can more easily emerge. 
20 This has been done using mechanical or ultrasonic cutting, chemical erosion (acidified 
Tyrodes solution), or by laser ablation of part of the ZP. 

A laser can be used to produce a trench in the edge of the ZP layer, penetrating 
through (or almost through) the ZP thickness. The trench creates a weakened region that 

25 provides a crack through which the embryo will later emerge. Lasers of many types can 
ablate the ZP. A laser of wavelength \ = 1480 nm, which is strongly absorbed in water, 
has been found to be effective for thermolysing the ZP. The laser can be used in pulses 
relatively short enough to avoid significant thermal conduction into the nearby embryo 
blastomeres, and at the sarne time avoid a chemical effect on the cellular chemistry, 

30 since it is in the non-actinic infrared wavelength region. The ZP is removed out to a 
radius around the laser beam determined by the local temperature history during the 
laser pulse. : 
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The same laser system can be used to ablate a.larger region of ZP so that an 
intact blastomere can be removed for external analysis. The laser is used, in this case, in 
a series of pulses directed at adjacent parts of the ZP to erode away a larger region. 
5 Typically, a gap is opened until a pipette can be introduced to suck out a blastomere. 
The embryo is relatively resilient, and generally survives both hatching and biopsy ZP 
ablation. 

Another related application of the laser system is in direct removal of the polar 
10 body for genetic analysis. Removal can be done either at the oocyte stage (first polar 
body), or after fertilization at the embryo stage (second polar body). Both polar bodies 
can be used to derive information on the genetic composition of the embryo. The 
procedure is analogous to laser assisted biopsy, except that in this case only the ZP layer 
and not the perivitelline membrane is penetrated, since the polar body remains between 
15 the perivitelline layer and the ZP. Additional related applications include transfer of part 
or all of the nucleus (nuclear transfer or transgenetic engineering), and ablation and 
destruction of part or all of the cell nucleus or oocyte spindle (e.g. in a cell to be used as 
the recipient of nuclear transfer). All of these applications benefit from the precise 
ablation capabilities of the laser system. 

The above applications are by way of example and should not be construed as 
limiting the possible uses of the invention. The present invention can be applied in any 
field where a laser beam is used with a microscope assembly. 

25 Summary of the Invention 

There is a need in the art for a configuration for introducing laser energy to a 
microscope, wherein the length of the laser beam emitted is reduced, and the quality is 
improved. There is also a need in the art for a method of calculating and displaying the 
30 isothermal contours associated with the energy generated in a sample by a laser. The 
present invention is directed toward further solutions to address these needs. 
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An optical injection system for use in conjunction with a microscope having a 

■ 

turret supporting an objective having an optical path is provided. The optical injection 
system includes a dichroic mirror disposed between the turret and the objective, along 
the optical path. A coUimating lens is further provided having an optical path directed to 
5 intersect the dichroic mirror. A laser source is positioned to project a laser beam 
through the coUimating .lens and along the coUimating lens optical path. 

In accordance with further embodiments of the present invention, the laser 
source includes an internal laser powered from an electronic package. The laser source 
10 can further originate using fiber optic cable providing optical energy from a laser. The 
laser can be one of an internal laser and an external laser, and operate at a wavelength 
(X) of about 1480nm. The laser beam can also be confocal with an optical image. 

The microscope can take the form of an inverted or a non-inverted niicroscppe. 
15 The microscope objective can take the form of a shortened optical train to provide a 
standard total length 45nMn parfpcal objective unit module. The system can also be 
adapted to longer parfocality systems, e.g., the Nikon CFI 60 optics (60 mm parfocality 
system). 

20 In accordance with another embodiment, the turret has a removable turret 

adapter. The turret adapter can be one of a plurality of turret adapter designs, thereby 
creating a universal mounting system, which enables the optical injection system to be 
mounted to the microscope. 

25 In accordance with further embodiments of the present invention, an electronic 

package can be disposed separate from the system. Alternatively, the electronic package 
can be built-in to a board in a computing apparatus. The laser can be mounted on a 
piston in a cylinder. Focus of the laser in such an arrangement can be provided by a 
single screw to adjust the distance from the laser to the coUimating lens. 

30 
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In accordance with another embodiment of the present invention, the laser may 
be located in the electronic package, and the laser energy transferred to the optical 
injection system through a fiber cable. 

The system of the present invention can be utilized in a number of different 
applications, including ablating, dissecting, moving, holding, or otherwise effecting 
biological cells or tissue with a laser source. 

A method for calculating and displaying the isothermal contours generated by a 
laser in a sample is also provided. The method includes applying a laser beam to the 
focal point of a sample, dividing the region near the focal point into cylinders coaxial 
with the beam, deriving the maximum temperature reached during the laser pulse of at 
least three points at arbitrary distances from the focal point, plotting the temperatures 
calculated as a function of distance from the focal point sufficient to generate isothennal 
contours, and generating a computer display of said isothermal contours corresponding 
to the temperature calculations. 

In accordance with further enibodiments of the present invention, the sample is 
placed in an isotropic medium and the temperature calculations are displayed as rings 
centered around the focal point. A picture of the sample can also be displayed with the 
rings. 

The method of the present invention can be utilized in conjunction with a number 
of different applications, including ablating, dissecting, moving, holding, or otherwise 
effecting biological cells or tissue with a laser source. 
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Brief Description of the Drawings 

The features and advantages of the present invention will become better 
understood with regard to the following description and accompanying drawings, 
5 wherein: 

FIG. lA is a diagrammatic illustration of a conventional fluorescent port laser; 
FIG. IB is a diagrammatic illustration of a conventional filter-cube laser; 
FIG. IC is a diagranmiatic illustration of a conventional filter-cube fiber optic 

laser; 

10 FIG. 2A is a side view illustration of a laser and microscope arrangement 

according to one aspect of the present invention; 

FIG. 2B is a side view illustration of the arrangement of FIG. 2A in further 
detail according to one aspect of the present invention; 

FIGS. 3A and 3B are perspective illustrations of a fiber optic and microscope 
15 arrangement according to one aspect of the present invention; 

FIG. 4 is an exploded diagrammatic illustration of a first embodiment of the 
laser assembly according to one aspect of the present invention; 

FIG- 5 is a flowchart illustrating the capture of images and conveyance to the 
computing apparatus according to one aspect of the present invention; 
20 FIG. 6 is a flowchart illustrating the generation of isothermal rings according to 

one aspect of the present invention; 

FIG. 7 shows the isothermal contours generated by applying a laser pulse of 
power 140 mW and duration 3 msec to the ZP of a pre-embryonic cell. The central ring 
shows corresponds to the laser Gaussian beam diameter, and the rings successively 
25 radiating from the central ring show the peak temperatures of 140, 100, 80, 60 and 50 
^'C, respectively. 

FIG. 8 shows laser ablation of ZP of bovine eggs at 100 mW, pulse duration 1.5, 
3, 4.5 ms. 

FIG. 9A shows bovine eggs in which several channels have been cut by laser. 
30 FIG. 9B sho\ys a close-up view of the channels in bovine ZP. The channels are 

approximately 25 jim long, almost constant-radius cylindrical intercepts with sharp 
edges and do not show effects of beam convergence or divergence. 
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FIG, lOA shows FEA analysis of 3 ms pulse, 100 mW, converging Gaussian 
attenuated beam solved in half plane. The isotherms at 0.1 ms into the pulse are 
indicated. Isotherm interval is 2.5^. Central temperature is 98.8^C. 

FIG. lOB shows FEA analysis of same case as in FIG. 9A, at 3 ms. Isotherm 
5 interval is 2.5°C. Central temperature is 

FIG. 11 shows GE analysis of flattop standard case (100 mW, 3 ms, a=3 |Lim), 
beam on axis of a 200 |xm right cylinder. Temperature given at various radial distances 
on focal plane. Using a Gaussian-profile beam increases the central temperature peak by 
about 15°C, but has a negligible effect on the temperature exterior to the beam. 
10 FIG. 12 shows the measurement of beam steering (horizontal view). 

FIG. 13 shows the bending of ray with impact factor b near a central refractive 
field (vertical view). 

FIG. 14 shows the probe beam steering angle in the IR laser thermal field versus 
power. Measurements are denoted by symbols and predictions by solid 
15 lines. The dashed line corresponds to the prediction when maximum temperature cannot 
exceed lOO^C. ; 

FIG. 15 shows the measured channel diameter versiis pulse duration averaged 
over three bovine eggs and two murine embryos. Predicted maximum temperature in 
focal plane at given radius is shown versus pulse duration and radial distance from axis, 
20 FIG. 16 shows the predicted maximum of the temperature pulse 20 \im from 

beam axis in the focal plane is given versus laser pulse duration and beam power. Also 
shown are liniit lines corresponding to beam central temperature 150 and 1 lO'^C. 

Detailed Description 

25 . 

An illustrative embodiment of the present invention relates to a system 
combining a laser with a microscope to result in an efficient and multi-functional device 
for performing micromanipulation of small objects with a laser. 

30 FIGS. 2A through 16 illustrate example embodiments of a laser assembly and/or 

a laser and microscope arrangement, in addition to example applications for use of the 
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laser and microscope arrangement, and results thereof, according to aspects of the 
present invention. Although the present invention will be described with reference to 
the example embodiments illustrated in the figures, it should be understood that many 
alternative forms can embody the present invention. One of ordinary skill in the art will 
5 additionally appreciate different ways to alter the parameters of the embodiments 

■ 

disclosed, such as the size, shape, or type of elements or materials, in a manner still in 
keeping with the spirit and scope of the present invention. 

Lasers have been utilized for manipulation of cells and cellular organelles. The 
10 • application of a precise dose of radiation for a set time in a pre-designated area can be 
used to destroy or neutralize organelles, such as the nucleus, as well as cutting, trapping 
and heating entire cells. The present disclosure utilizes an example of using the lasers 
for drilling into embryo zona pellucida ("ZP"). However, this is merely an illustrative 
application and in no way is it intended that the present invention be limited to such 
15 application. 

T 

Adapters and kits used in conjunction with microscopes are shown in two patents 
described below. US Patent No. 5,349,468 is directed to an adapter for a microscope 
that fits between the objective and the microscope turret for providing fluorescence 
20 microscopy. In this device a light source provides illumination in a direction normal to 
the optic axis, which projects the light after collimation onto a dichroic mirror, and then 
to the objective, to give epifluorescent illumination of the target. 

US Patent No. 4,884,880 is directed to a "kit" that converts the standard 
25 microscope into a single aperture confocal scanning epi-illumination microscope. The 
kit fits between the objective and the microscope for the purpose of creating scanning 
epi-illumination. 

In a conventional fluorescent port microscope arrangement, laser radiation is 
30 conveyed from the source to the microscope objective using a collimated beam and an 
on-axis dichroic mirror (see FIG* lA). The optical system is analogous to an 
epifluorescent design, in which laser radiation is used in place of the exciting fluorescent 
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illumination. Chromatic aberration between the laser infrared (IR) beam and visible 
light is removed by adjusting the degree of coUimation of the laser beam, so that it is 
confocal with the visible image and produces its maximum effect exactly in the same 
plane as the image. In some conventional devices, long IR pathlengths were used to 
5 convey the laser energy between a laser source 10, generally at the back of the 

microscope, through a fluorescent cavity 12, on to a dichroic mirror 14, and through an 
objective 16 to the target. The objective 16 must be highly transnaissive in both the laser 
radiation wavelength and the visible wavelength for such an arrangement. The long path 
length resulted in attenuation of the laser beam and reduced its safety, since prevention 
10 of local blastomere heating is more difficult with less powerful beams. A typical 
arrangement of this type is shown in FIG. 1 A. 

More recent conventional designs, as shown in FIG. IB, minimize the 
pathlength between a laser source 20 and an objective 22 by packaging an entire optical 

15 system 24, including a dichroic mirror 26, into a region under the objective 22. This is 
the same location where the ifluorescent filter-cubes are placed in fluorescent operation. 
This arrangement increases the laser power available and makes the laser more portable 
and easily adjustable. However, the arrangement prevents use of fluorescent 
illumination, or other optical devices needing the same space (e.g. the Spindle View^M 

20 polarizer, made by GRI, Inc. of Wobum, MA, for detecting the oocyte spindle for ICSJ), 
and the laser package generally must be removed before devices of such type can be 
utilized. 

In certain other conventional designs, the IR radiation is conveyed to an 
25 objective 30 using a fiber-optic cable 32 as shown in FIG. IC. The exit of the fiber- 
optic cable 32 acts as a laser source, and the beam passes through a coUimating lens 34, 
which may be, for example, an aspheric or a GRIN lens, and reflects from a dichroic 
miitor 36 to provide the required almost-parallel beam into the objective 30. This 
arrangement has the advantage of facilitating the positioning of the fiber-optic "source". 
30 Moreover, another (e.g. a visible) laser cable can be spliced into the fiber-optic cable 32, 
giving in effect tv^o coincident sources. However, the transfer of energy into and out of 
the fiber optic cable inevitably reduces the laser power available, due^to focusing. 
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reflective, and cable attenuation losses. The available energy therefore tends to be lower 
in fiber optic systems, again leading to the use of longer pulses and the associated 
greater heating of neighboring blastomeres. 

One difficulty associated with laser power delivery by fiber optic cables is that, 
due to the very small diameter of the fiber, the presence of tiny particles or other 
obscuration at the cable entry or exit can greatly diminish the intensity of the laser beam 
dehvered to the target. While this is a common issue for the conventional digital 

¥ 

information transfer application of fiber optic cables, such digital signals can sustain a 
loss and still transmit information properly. However, it is a greater potential difficulty 
with energy-transfer applications where the unreliability of the beam intensity at target 
makes accurate and repeatable ZP hole-drilling problematic. 

The disadvantages of the prior methods can be overcome by the teachings of the 
present invention. In accordance with the present invention, advantage is taken of the 
fact that a laser source is packaged in a spiall integrated-circuit "can" 66 shown in FIG. 
2B, and is sufficiently small that it may be positioned very close to an objective 54. In 
one embodiment of the present invention, a laser assembly 76 including a specially 
designed low-profile objective (e.g. a short. 21 mm objective) is combined with an 
optical injection system 56 to introduce the laser light between a microscope turret 58 
and the objective 54. The turret 58 (see FIG. 2A) includes a removable turret adapter 59. 
The turret adapter 59 is one of a number of different turret adapter designs. The 
plurality of turret adapter designs creates a universal mounting system, which enables 
the laser assembly 76 to be mounted to any microscope. 

In this design, a fluorescent filter-cube channel 60 (see FIG. 2A) beneath the 
inverted niicroscope turret 58 is left unoccupied, and may be filled with any optical 
equipment required, since it is no longer integral to the optical design of the laser 
system. The entire laser assembly 76 exists between the microscope turret 58 and a 
microscope stage 62 as shown in FIG. 2A. The position of the laser assembly 76 
between turret 58 and stage 62 is indicated. 



HTS-035DV 

11 

The optical injection system 56 is shown in additional detail in FIG. 2B. A 
solid-state laser beam is directed through an adjustable coUimating lens 68 on to a 
dichroic mirror 70. The dichroic mirror 70 reflects the wavelength X = 1480 nm 
radiation beam along the optic axis and up into the objective 54, in a direction opposite 
5 to the image beam light. A laser diode 66 is mounted inside a piston 72 of precise 
dimensions, which slides in a cylinder 74. The sliding movement of the piston 72 
controls the distance of the laser diode 66 to the coUimating lens 68, while preventing 
lateral movement. Any chromatic aberration between the laser beam and the visible 
image is removed by adjusting the piston 72 position. The movement of the piston 72 
10 focuses the laser beam. 

The laser diode 66, piston 72, coUimating lens 68 and dichroic mirror 70 are 
packaged in optical injection system 56. The objective 54 screws in to the top of the 
optical injection system 56, creating the laser assembly 76, which can be coupled to the 

15 turret 58 by way of the turret adapter 59. The entire laser assembly 76 is constructed so 
that it can be mounted like an objective, by first screwing the tuijet adapter 59 on to the 
objective nosepiece turret 58 and then screwing the laser assembly 76 on to the turret 
adapter 59. The design of the laser assembly 76 to turret adapter 59 interface is such 
that the laser assembly 76 may be oriented in any rotational position and then locked 

20 down into place. This feature overcomes the problem of a fixed thread design, which 
would have only one final rotational position. In the example embodiment, the total 
length of the laser assembly 76 including the length of the specially designed objective 
54 is equal to the length of a standard 45mm parfocal objective. The complete laser 
assembly 76 fits between the turret 58 and the stage 62 of a microscope. Although the 

25 system is designed for use in inverted microscopes, it can easily be used in upright 
microscopes as well. 

In the particular application of ZP hatching and embryo biopsy, a shortened 
objective 54 is utilized, with magnification 50x and NA = 0.6. The objective 54 is made 
30 to transmit more than 90% power at the wavelength of A, = 1480 nm. The laser diode 66 
is precisely positioned within the laser assembly 76, and is maintained in the required 
rigid optical relationship with the other parts of the laser assembly 76. The optical 
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injection system 56 therefore remains with the required focal and directional alignment 
between visible and infrared beams. The entire laser assembly 76 can therefore be 
transferred to other microscopes 65 {see FIG. 2A) without changing the relationship 
between the optical components forming the laser assembly 76, and without requiring 
realignment. 

Electrical power to the laser diode 66 is provided from a standard control circuit 
78 on a conventional printed circuit board (PCB) 80. The PCB 80 may be mounted in a 
separate box, or as a board inside a computing apparatus 82. A cable 84 between the 
computing apparatus 82 and the laser assembly 76 connects each component to the 
other. The cable 84 normally connects to the laser assembly 76 with a rotational swivel, 
such that the microscope turret 58 {see FIG. 2A), and any of the other objectives used, 
can rotate without causing inconvenience from the controlling cable 84. It should jbe 
noted that inclusion of Zener and Schottky protection diodes (not shown) mounted on 
interface PCB 106 of FIG. 4 within the laser assembly 76 minimizes risk of standard 
electrostatic-damage to the laser diode 66. The laser diode 66 is further protected by 
automatically shunting the laser diode 66 connections when the control cable 84 
disconnects from the interface PCB 106 within the laser assembly 76. 

A software application installed on the computing apparatus 82 controls the laser 
assembly 76 (laser power, laser pulse length, laser multipulse). The phrase "computing 
apparatus" as used herein refers to a progranmiable device that responds to a specific set 
of instructions in a well-defined manner and can execute a set of instructions. The 
computing apparatus 82 can include one or more of a storage device, which enables the 
computing apparatus to store, at least temporarily, data, information, and programs (e.g., 
RAM or ROM); a mass storage device for substantially permanently storing data, 
information, and programs (e.g., disk drive or tape drive); an input device through which 
data and instructions enter the computing apparatus (e.g., keyboard, mouse, foot switch, 
or stylus); an output device to display or produce results of computing actions (e.g., 
display screen, printer, or infrared, serial, or digital port); and a central processing unit 
including a processor for executing the specific set of instructions. The computing 
apparatus has access to code for transmission of the niicroscope image on to a monitor 
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(not shown). The code includes routines for measuring the image, storing and retrieving 
electronic images, operating the laser assembly 76, and generating and printing reports, 

s 

As shown in FIGS. 3A and 3B a fiber-optic cable 86 may be substituted for the 
laser diode 66 in an optical injection system 57. A laser 77 is then located at the other 
end of the fiber optic cable 86, together with an electronic laser control circuit 79. This 
may be built in to an optical injection system 57, in a separate system, or built on to a 
PCB for insertion into a computing apparatus. In the latter case, the fiber-optic cable 86 
plugs directly into the PCB in the computing apparatus and no electric connections are 
required between the assembly and the computing apparatus. In the example 
embodiment, the laser 77 provides the necessary signal through the fiber-optic cable 86 
to the dichroic mirror 70, which reflects the signal to the objective 54 as before. The 
substantial difference in this embodiment is the use of the fiber optic cable 86, which 
can allow the laser 77 source to be located in a number of locations proximal or distal 
from the optical injection system 57. 

The use of fibers has the advantage of producing an almost Gaussian beam, 
allowing relatively good focusing properties. Fibers may be spliced to introduce a 
visible beam (for aim spot or focus spot) co-directional with the JR beam. However, 
fibers also have the same problem of particle obscuration mentioned above. 

It should be emphasized that the design used in the present invention can be 
applied to substantially any embodiment of a laser to be used for microscope irradiation. 
Some examples include laser cutting, e.g., at a wavelength of A, = 337 -390 nm, laser 
scissors and laser tweezers, and laser differential heating. DNA denaturing (e.g. nucleus 
ablation) is another application. The system can be extended to any laser from which 

the power is delivered in a fiber-optic cable. The present invention is not limited to 

> ■- ' > 

inverted microscopes. It can be used in an upright microscope, and requires no , 
adjustnients when transferred between one type of microscope and another. 

FIG. 4 illustrates another detailed example embodiment of an optical injection 
system 88 built in accordance with the teachings of the present invention. A first 
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embodiment of the optical injection system 88 was designed and machined, and is 
shown in FIG. 4. A body 90 houses the elements of the optical injection system 88. A 
mirror mount 92 provides support for a dichroic mirror 102. An adapter nut 94 attached 
to the body 90 with a retainer ring 114 mounts the laser system to a turret adapter 96. A 
diode mount 98 and diode clamp 100 support a laser diode 120, while a compression 
spring 1 16 pushes against the diode mount 98. A focus screw 104 is used to adjust the 
laser light. A connector on interface PCB 106 is suitable for connecting the assembly 
with a computing apparatus (such as the computing apparatus 82 of FIG. 2B). A lens 
bushing 108 positions a lens 118. A cover 1 10 surrounds the laser portion of the optical 
injection system 88. A compact objective 1 12 mounts substantially orthogonal to the 
laser diode 120 direction. 

The optical injection system 88 was used with a specially designed HTB short 
lens, and the assembled unit with lens had a parfocallength of 45 mm, equivalent to the 
length of a standard objective. The optical injection system 88 was demonstrated on a 
Nikon TE 300 microscope, then transferred to a Leica DMEL, Zeiss Axiovert 25 and an 
Olympus IX-70 microscope. In all cases, the focus of the IR beam did not have to be 
adjusted, and the unit gave high quality images. Standard tests on the HTR B^im grid 
showed no detectable image difference between the microscopes. 

r I 

The IR transmission is also independent of the microscope, since it can only 
depend on the optical injection system 88. The value of the transmission, T, was 
measured with the system by removal and replacement of the objective, receiving the 
beam in an OPHIR Nova NIST -traceable detector. The result was T = 90%, close to the 
values obtained with the standard objective. 

The size and quality of the holes drilled in test bovine oocyte samples was 
measured. The results fell within the predictions expected from the thermal diffusion 
transfer theory [Douglas-Hamilton DH and Conia J (2001)]. This result would be 
expected since the entire imaging apparatus is contained within the optical injection 
system 88. 
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The teachings of the present invention provide for JR focal adjustment to be 
required only at installation. The laser assenably can be fitted to any microscope without 
changes. The optical injection system works on inverted and on upright direct 
transmission microscopes. Short optical path gives maximuni beam power available. 
5 High beam power avoids damaging local cells that are not targeted. The optical 

injection system can be adapted to fiber optic sources. The entire electronics associated 
with the optical injection system can be contained inside computing apparatus. The 
optical injection system is small and lightweight, and does not require special 
installation. The optical injection system can further be returned for factory servicing, 
10 and re-installed by a user without losing its alignment or requiring optical adjustments. 

An illustrative embodiment of the present invention also relates to a method for 
calculating and displaying the isothermal contours generated by a laser in a sample is 
also provided. As used herein, the term "isothermal contours" refers to the array of 

15 maximum temperatures generated by a laser pulse at particular locations within a 
sample, e.g., the ZP of a pre-embryonic cell. The locus of all points reaching the 
maximum temperature represents an isothermal contour of maximum temperature or 
peak thermal excursion. The method includes applying a laser beam to the focal point of 
a sample, dividing the region near the focal point into cylinders coaxial with the beam, 

20 deriving the maximum temperature reached during the laser pulse of at least three points 
at arbitrary distances from the focal point, plotting the temperatures calculated as a 
function of distance from the focal point sufficient to generate isothermal contours, and 
generating a computer display of said isothermal contours corresponding to the 
temperature calculations. 

25 

The isothernial contours illustrate the thermal effect of a laser beam on a sample, 
e.g., the ZP of a pre-embryonic cell. A wavelength of approximately 1480 nm radiation 
is strongly absorbed in water and, therefore, absorbed in all living tissue. Consequently, 
the region of the laser beam in an aqueous medium becomes significantly heated and the 
30 temperature rise at the focal point can easily reach AT = 150 °C by the end of the laser 
pulse. The ZP melts at a temperature of about T = 140 ''C. Therefore, the ZP will be 
removed within a ring corresponding to peak temperature T of about 140 °C. 
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The thermal energy deposited by the laser beam is conducted away through the 
medium. During the laser pulse a dynamic balance is set up between the laser heating 
and the thermal conduction rate which results in the temperature rising rapidly until a 
5 steady-state temperature is reached. Generally, the pulse is terminated long before the 
steady-state. Consequently, regions at different distances from the focal point reach 
different temperatures that depend on the pulse duration, the beam intensity, distance 
froni focal point, and conductivity of the medium. The temperature at or near the focal 
point of the beam can be as high as AT > 200 "^C. The temperature falls rapidly with 
10 increasing distance from the focal point. In order to ensure the safety of the laser system 
and the success of the laser beam manipulations, it is very important to accurately 
calculate the particular temperatures at various positions of the sample. 

The medium containing the sample should be compatible for various types of 
15 organic structures, such as embryos, blastomere, or other cells. These structures will 
differ in chemical composition from the medium itself and in general will have different 
thermal properties. The thermal diffusion equation can be solved in this type of 
geometry and the thermal history derived at every point. The locus of all points 
reaching a particular maximum temperature represents an isothermal contour of peak 
20 thermal excursion. However, since the major component of most biological structures is 
water, the samples will have thermal and IR-absorption properties similar to those of 
water. 



A theoretical derivation of the thermal diffusion regime for an isotropic medium 
25 has been given by Douglas-Hamilton and Gonia (2001) (/. Biomed Opt (2001) 6:205-13 
"Thermal effects in laser-assisted pre-embryo zona drilling")* incorporated herein by 
reference. The region near the focal point is divided into cylinders coaxial with the 
beam and the thermal history at arbitrary distances from the focus is derived which 
determines the thermal regime or isothermal contour as a function of distance. Such 
30 values also specify the distance from the focus required to maintain safety. In particular, 
these calculations are applied to cells that are within the embryo while the ZP is being 
drilled, in order to maintain a safe stand-off distance between the drilling operation and 
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the embryonic cells. The temperature of the medium near the beam was measured 
experimentally and confirms the theoretical thermal predictions. Any mediuni 
anisotropy can be included in the calculation and the appropriate isothermal contours 
derived. 

5 

Using the foregoing thermal profile, the maximum temperature reached during 
the laser pulse as a function of distance from the beam axis is plotted and expressed as a 
series of isotherms, each representing the peak temperature reached during (or after) the 
laser pulse. To first approximation, the medium may be treated as isotropic and the 

10 peak-temperature isotherms contours are generated as rings centered on the beam axis or 
focal point. The isotherms are displayed on a computer display centered on the beam 
focal point and can be identified by color or label to give the local peak temperature 
computed from the thermal diffusion model. Because the local heat input depends on 
the laser power and the laser pulse duration, the diameter of the isothermal rings will 

15 vary as a function of the laser power and pulse duration. 

The computing apparatus 82 can display the isotherms as described in the above 
example application, or a different computing apparatus can be utilized to display the 
isothenris. Prior to displaying the isotherms, the computing apparatus 82 must first 

20 obtain images of the sample as viewed through the microscope. To accomplish this, and 
as illustrated in FIG. 5, the microscope images are either taken with a digital camera or. 
an analog CCD camera (step 200). If the images are taken with a digital camera, or are 
otherwise in digital format, the images are transferred directly onto the monitor screen 
(step 202). If the images are not in digital format, the analog signal is transferred to a 

25 digitizer board (step 204) to be digitized [e.g. the Bandit digitizer/display board from 
Coreco Corp, Quebec, or equivalent]. The digitized image is then transferred directly to 
the screen (206). In either instance, the image can be over- written with information 
from a digitizing system, so that any computer information may be superimposed on the 
image. 

30 

Once the image of the sample is captured from the microscope, the computing 
apparatus 82 takes the image and superimposes a plurality of isotherms in the form of 
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lines drawn on the image of the sample. The process carried out by the computing 
apparatus 82 for creating the isotherms can be described as follows, and as illustrated in 
FIG. 6. 

r ' ' 

■* .-r * 

5 In the example embodiment, the Isotherms are generated by interpolating values 

from fixed solutions for a 180 mWatt beain of radius 2 micron. The fixed solutions are 
computed at pulse durations of 0.5, 1,2,3, 4, 5, 10, 20 milliseconds (step 210). The 
isotherms for a beam of power 180 mW at any given pulse duration can be accurately 
interpolated from the fixed solutions (step 212). By way of example, using a third 
10 degree polynomial interpolation gives results within 5% of the correct values: Since the 
temperature excursion scales directly with beam power, only the coinputed solutions for 
different pulse durations at a particular given beam power are required, The 

■ * 

corresponding isotherms can then be computed at any other beam power, for the same 
beam radius, by direct linear extrapolation (step 214). 

15 ■ ^ 

In practice the system, is used at maximum power, and the effect of the beam is 
controlled by decreasing the laser pulse duration. In this way, the conduction of heat 
from the laser beam to local (e.g. blastomere) cells is minimized. 

20 The corresponding peak temperature isothermal ring diameters are derived from 

the laser beam power and the laser pulse duration input at the controls (step 216). For 
any cornbination of pulse duration and beam power, there is a single-valued set of rings, 
corresponding to the positions at which each peak temperature is reached. 

25 It is of great advantage to the user of a laser manipulation system to be able to 

see at a glance what the effect of the laser will be at positions remote from the actual 
laser focal point. Therefore, the derived isothermal contours are displayed on a 
computer screen. In an isotropic medium the contours are rings. The ring diameter 
depends on the laser power and the laser duration. As the laser power and the laser 

30 duration are changed by the user, the rings expand or contract accordingly. This enables 
the user to see exactly what the potential for embryo damage is with the laser, and to 
adjust the parameters accordingly. 
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In application to the drilling of ZP, it is a further advantage of the invention that 
the size of the hole to be drilled is represented well by the diameter of the 140 ""C peak- 
temperature isotherm. Thus, for any conditions of laser power and laser pulse duration, 
5 the diameter of the Tmax = 140 ''C isotherm gives the diameter of the hole that will be 
drilled in the ZP. 

An example of the invention is shown in the FIG, 7, in which the isotherms are 
shown superimposed on a scale image of an oocyte indicating the peak thermal 
10 excursion that will be experienced at different distances from the laser beam focus by 
different parts of the oocyte. 

r 

1 

The present invention is further illustrated by the following examples which 
should not be construed as further limiting. Numerous modifications and alternative 

15 embodiments of the present invention will be apparent to those skilled in the art in view 
of the foregoing description. Accordingly, this description is to be construed as 
illustrative only and is for the purpose of teaching those skilled in the art the best mode 
for carrying out the present invention. Details of the structure may vary substantially 
without departing from the spirit of the invention, and exclusive use of all modifications 

20 that come within the scope of the appended claims is reserved. Such equivalents are 
intended to be encompassed by the following Examples and claims. 
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Examples 

Example 1 : Laser Module, Microscope Specifications, and Design 

A Nikon TE-300 laser-equipped inverted microscope was used, which has the IR 
5 laser beneath the stage. Laser light proceeds up through the objective and is focused on 
the target. The laser module (ZLTS, Hamilton Thome Research, Beverly, MA) has 
maximum power 200 mW at k- 1480nm, and is utilized in pulsed mode. Peak 
transnMtted laser power through the objective is approximately 140 mW. Pulse duration 
is adjustable in half millisecond increments from 0.5 to 25 ms. In addition to the laser 
10 diode itself, the laser module includes control board, adjustable coUimating lens and 
dichroic mirror, and is inter- locked to prevent potentially hazardous operation. 

The laser module fits into the microscope filter cube slot under the stage and the 
beam is deflected by a dichroic mirror along the optic axis of the microscope. The 

15 nearly coUimated beam is directed toward the back aperture of a 40X, 0.60 numerical 
aperture (NA) objective lens, designed to maximize transmission at 1480 nm. The long 
working distance (WD = 3 mm) objective lens is used to focus the laser beam on to the 
target specimen; it has high transmission ratio in the infrared range (typically about 71% 
at 1480 nm). The collimation of the beam is adjusted to give parfocality of DR. and 

20 visible light. Biological specimens, e.g., bovine oocytes or mouse pre-embryos, are 
maintained in an aqueous culture medium in clear polystyrene Petrie culture dishes 
(Falcon, Lincoln Park, NJ). The incident focused laser beam travels sequentially 
through air, plastic, and then the aqueous solution before reaching the target area of the 
zona pellucida (ZP). 

25 

Delivered laser energy at the objective was measured using a thermal detector 
head (Ophir Optronics, Peabody, MA) placed immediately above the stage with 
objective removed. Transmission of the objective was determined by placing an 
identical opposed objective above the focal point to recoUimate the beam which passed 
30 into the detector. Using the derived transmission, 71 % at X = 1480 nm, maximum beam 
power at the focal point was estimated at 103 mW, allowing for reflection at the 
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polystyrene Petrie dish and absorption through 75 iim H2O. Pulse duration was checked 
using an oscilloscope (Tektronix TDS 360) monitoring laser radiation scattered on to a 
Ge photodiode detector (B 1918-01, Hammamatsu, Bridgewater, NJ). In a preliminary 
measurement, the focal beam radius in air was determined by intercepting the beam at 
the focal point with a knife edge which was moved across the beam at known rate while 
the IR intensity was monitored with the Ge detector. At the focus, the radius to the 1/e^ 
intensity level was R=3.1 ± 0.5 jim. The focal cone half angle in the medium is 26.7^. 

Example 2: Laser-Assisted Zona Dissection v^^ith Bovine Eggs 

Bovine eggs are readily available. They constitute a useful and cost effective 
model for testing laser-assisted zona drilling protocols. The Nikon TE300 inverted 
microscope is used with illumination source and condenser above specimen and 
objective below. The Petrie dish specimen chamber (wall thickness 1 mm, diameter 50 
mm) is on a motorized stage. The eggs settle at the bottom of the dish and the objective 
is focused on the edge of the ZP, about 75 fim into the medium. 

1 

A typical zona drilling is shown in FIG. 8 for a bovine oocyte. The egg diameter 
is about 150 ^in and the thickness of the zona pellucida, the translucent shell 
surrounding the egg, is about 12 fim. The laser power was set to 100 mW delivered at 
the focal point and three points in the ZP corresponding to pulse duration 1.5, 3, and 4.5 
ms. Beam diameter to 1/e^ was 6.2 \im, independent of beam power. The hole 
diameters are 4.3, 7.8, and 10.5 \im; the hole diameter depends more on pulse duration 
than on laser beam size. 

In FIG. 9A an electron photomicrograph of several channels in a bovine egg is 
shown. Repeated firing at constant laser pulse length 25 ms and laser power 50 mW 
produced channels of consistent size. 

The close-up in FIG. 9B emphasizes the sharply drilled wall of each channel. 
The channels have clear cylindrical cross sections of width approximately 12 jxm with 

no significant beam broadening at either end. The channels reflect a thermal dissolution 

/ 
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"melting" temperature and its locus is. not coterminous with the beam edge. Although no 
evidence of thermal damage can be seen surrounding the target area, the transient 
thermal excursion may cause damage invisible under the microscope to nearby living 
cells. The thermal history of the laser drilling is derived as described below. 



10 
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20 
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EXAMPLES: Thermal Predictions 

Standard Conditions 

The objective (NA=0.6, 40x, infinity corrected) used transmits 71% of the 
incident energy at X = 1480 nm. The typical laser beam power arriving at the vicinity of 
the ZP following transmissive, reflective, and attenuation loss is calculated as 100 mW. 
The typical beam focal radius is calculated as a = 3 p.m, assuming the beam is passing 
up through an inverted microscope into medium of refractive index n= 1.333. The 
thermal histpry or isothermal contour expected along the path of the laser beam and 
especially in the vicinity of the pre-embryo is derived as described below. 



Thermal Constants and Source Function Geometry 

Using medium comprised of almost pure water and the cellular material is 
approximately 80% water, the infrared absorptances in the medium and cell at X = 1480 
nrn are approximately that of H2O and cannot be significaiitly lower. Similarly, the 
thermal conductivity of the material will not be higher than that of water. By taking the 
absorptance as identical to that of H2O, a = 21 cm"\ a lower limit to the temperature 
produced is obtained. On the time scales of interest (1-25 ms), the only significant mode 
of heat loss from the heated liquid will be thermal conduction, since convection does not 
have time to develop and radiation may be neglected. The thermal conductivity of water 
is 6 X 10" w/cm/° at room temperature, but since the thermal excursion during the laser 
pulse is typically 100 - 200''C, we take the conductivity of the medium as K = 6.8 x 10'^, 
corresponding to H2O at lOO'^C. 



The laser path is approximated as three regions: (1) a con verging cone of 
30 semiangle 0; (2) the waist region of the beam; and (3) a diverging cone with semiangle 
0. The beam has higher net power at the converging than at the diverging cone, due to 
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attenuation. The cone semiangle is 0 = arcsin (NA/n) = 26.7^ . For the thermal 
calculation we approximate the waist region as a cylinder of radius a = 3 |xm and length 
2a • cot 9. 

The laser beam converges through the medium to a focal point 75 ^im above the 
floor of the Petrie dish. Absorption will only be important in the medium and may be 
neglected in the Petrie dish. Self-focusing of the IR beam will not significantly affect 
the thermal distribution outside the beam and is ignored. The radiation intensity is 
symmetric in azimuthal coordinate, varying in intensity only with axial distance and 
radius. Ignoring angular variation, the heat diffusion equation in axisymmetric 
cylindrical coordinates may be written 

Equation 1: 

dT a^T I dT d^T S 
dt dr- r dr dv pC^ 

where the heat diffusion coefficient is k = K/pCp, and pCp = 4.18 j/cm^, with p the liquid 

i 

density and Cp its specific heat, so k= 1.6 x 10"^ cm^/s. Variation of the diffusion 

r 

coefficient with temperature is small and is ignored. The origin is taken as the focal 
point, with r and z as the radial and axial distances. The source function S represents the 
laser heating power per unit volume and varies across the entire domain. . 

■ 

The source functions in the three regions of the laser path are normalized to the 
laser power P reaching the focal plane so that the converging beam is more intense 
below the waist and attenuates as it progresses upward, diverging above the waist. The 
central cylinder length is short (2a.cot 9=12 |jim), and attenuation is ignored in region 2. 

The source functions for attenuated uniform beam in the three regions are 

■ » 

Equation 2 - Converging: 
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Equation 3 - Waist 
aP 

S2=- 



Equation 4 - Diverging; 



P3 = ' ---^ r.^ 

27i{l-cos0) r^^z^ 



10 The Gaussian profile is approximated by applying the factor a^"^^^^'^^^"'^^ ^^^^^ to Si 
and S3, where A is a normalization constant, and for S2, we use the factor Ae~^^^^^^^ . 



Boundaries for Solution of the Heat Diffusion Equation 

Equation (1) was solved by finite element analysis (FEA) using the codes 
15 COSMOS (Structural Research Co., Los Angeles, CA) and Mathlab PDE Toolbox 

(Mathworks Inc., Natick, MA) with source functions from Equations (2)-(4). The beam 
travels parallel to the axis of a cylinder of radius 100 [im and length 200 [im. For short 
pulses (x < 6 ms) these dimensions are sufficient to make edge effects negligible since 
the thermal diffusion distance in time x is x < 30 |Lim. The axis is a Neumann boundary 
20 with zero radial temperature gradient. Since the beam is attenuated, the focal plane 

cannot be taken as a Neumann boundary. Therefore, the equation must be solved in the 
entire upper half plane. The cylinder end planes at z = ± 100 |Lim and cylinder wall at r = 
100 \im are both Dirichlet boundaries with temperature held at To = SV'C. The initial 
temperature of the entire domain is taken as To. 



Of practical interest is thermal history on the focal plane itself at various radial 
distances from the axis, which will give information on the thermal excursion 
experienced by the nearest blastomeres. 
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Preliminary VERincATiON 

\ 

For comparison with an analytic case, the FEA was run with a uniformly heated 
cylinder of radius 3 Jim in a 200-|Jim -diameter right cylinder containing H2O at To = 37, 
and with source function as in Equation (2), P= 100 mW and a = 21 cm \ 



The analytic steady-state solution for the temperature difference between edge 
and center of an infinite cylindrical region of radius Rmax heated on axis by a uniform 
beam of radius R is 



10 Equation 5: 



dP 



AnK 



l + 21n( 



R 



max 



J 



The FEA steady-state result was within 0.5% of the temperature difference from 
15 Equation (5). Therefore, the systeni is integrating correctly. 



Standard Case 

The initial temperature is taken in the medium at the physiological To = 37 °C in 
a cylinder of diameter 200 ^im and length 200 jim, with attenuated converging beam of 

20 NA=0.6. The field is divided into 1.13 x 10^ elements. The same temperature To is held 
at the end planes and at the radius Rmax = lOQiim. Beam power is 100 mW at the focal 
waist, with pulse duration 3 ms and radius 3 [Jim. The problem is solved in the half 
plane and has been converted to the full plane in FIGS. lOA and lOB where the beam 
direction is bottom to top, shown for times 0.1 and 3 ms, respectively. At 0.1 ms the 

25 central temperature is 98.8°C. In FIG. lOB, at 3 ms the central temperature is 189°C. If ^ 
the beam is flattop instead of Gaussian, central temperature is reduced approximately 
15%. In either case, these high temperatures imply the potential presence of superheated 
medium. The temperature change from the initial value at anytime depends on the 
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pulse duration and is directly proportional to the beam power (see Equations (l)-(5)). 
Predictions for any beam power may be scaled from the results reported. 

The effect of attenuation is apparent in the slightly non vertical slope of the 
5 isotherms as the beam progresses from below to above the waist in FIGS. lOA and lOB. 
Apart from this effect, the converging and diverging parts of the beam away from the 
end planes result in almost cylindrical isotherms. For the present purpose the focal 
plane is the only region of interest in questions of cell heating. Therefore, the question 
of whether the focal plane temperature heating can be approximated by a simpler (and 

10 faster) oncrdimensional solution is examined. 

•I ' . . ■ 

To test the effect of geometry the FEA system was run for the Gaussian beam in 
the following three cases, all with power 100 mW, pulse 3 ms, and radius 3 |Jim: 

1. Attenuated beani, NA=0.6 (converging cone, central cylinder radius 3 nm, 
15 converging/diverging cone); 

2. Unattenuated beam, NA=0.6 (same as above, exponential terms in Equations 
(2) and (4) set to unity); and 

3. Unattenuated beam, NA=0 (uniform central cylinder radius 3 |bim, 
exponential terms set to unity). 

20 

The latter case corresponds to a constant cylindrical beam and has one- 
dimensional symmetry, whereas the first two require two dimensions. 

For all three cases, predicted temperatures on the focal plane for radii 5 < radius 
25 < 50 iim as function of time are within 0.5°C over 10 ms integration time. At least on 
the focal plane, in the present case the converging and diverging beam produces heating 
very similar to that from a long heated cylinder, and if the power is normalized to its 
value at the focal plane, the effect of beam attenuation may be ignored. Hence, the 
niathematically simpler case of the cylinder can be used for predicting the temperature 
30 history at points on the focal plane. Accordingly, the one-dimensional Gaussian 
elimination solution has been used as set forth below. 
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Example 4: Solution in Cylindmcal Symmetry 

, Gaussian Elimination 

The FEA solution shows that radial temperature history is almost cylindrical 
with isotherms close to parallel to the axis over the region of interest. If the heating 
effect of the beam is approximated as a cylindrical region of constant radius "a" equal to 
its focal radius, Equation (1) can be written as a one- dimensional difference equation, 
which can be rapidly solved by Gaussian elimination ' to give the thermal time 
history in the region near the laser focus. The advantage of this is that solution is much 
faster than FEA, and the case of pure cylindrical geometry without converging and 
diverging beams gives a good approximation to the focal plane temperature. A Gaussian 
eUmination (GE) code to model the thermal behavior in the radial dimension is 
described herein. 

The area 0 < r < Rmax is covered and the Rmax = 100 |Lim is chosen because at the 
laser pulse lengths considered, the thermal excursion does not reach Rmax before the end 
of the pulse. The cylindrical area is divided into N shells, each separated by Ar = 
Rmax/N. Accurate results are obtained with Ar = 100 |Lim. As boundary conditions, the 
initial temperature and the temperature at r = Rmax are set to To = 37°C. The temperature 
gradient at beam center is set to zero. In the GE analysis, it is assumed that the beam is 
flattop rather than Gaussian: this reduces central temperature by about 15%, but does not 
affect temperature history outside the beam. Equation (5) is used to check the result and 
find agreement within 0.1% between the GE code and Equation (5) for very long (200 
ms) pulses, so the GE code gives correct results for the analytic case. 

Thermal Histories on Focal Plane 

The GE solutions for temperature at various radial positions on the focal plane 
are shown as function of time in FIG. 11. The standard case of beam power 100 mW, 
pulse duration 3 ms, in a 200 |Lim right cylinder initially at 37°C is used. At beam center 
the peak temperature reaches over 170°C. Temperature falls off sharply as the pulse 
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terminates due to the short diffusion time over 3 Jim. As expected, the thermal 
excursion seen near the beam has approximately the same duration as the laser pulse. 

Superheating , 

5 The calculated (GE) peak central temperatures at typical beam powers with beam 

radius 3 ^im in water are given for flattop beam focal spots of various pulse durations in 
Table 1 with steady-state values from Equation (5) for reference. Again, the region of 
the analysis is taken as a cylinder of length 200 ^im and radius 100 |Xm. The calculated 
temperature of water exposed to the beam is very high, with peak central beam 

10 temperature ranges up to above 200°C, corresponding to highly super-heated water. The 
question is whether this temperature is real or whether a phase change occurs which 
would reduce the central temperature and form a column of vapor bubbles along the 
beam axis near the focal point. Miotello and Kelly^^ have examined the formation rate 
of homogeneous nuclei and subsequent explosive phase change in superheated liquids. 

15 In the absence of a surface for heterogeneous nucleation, homogeneous nucleation is 
necessary for rapid phase change, and its rate depends strongly on how close the liquid 
is to the critical temperature (Tc = 647 K for H2O). The homogeneous nucleation rate 
becomes significant near 300°C. However, at and below 250°C, the rate of 
homogeneous nuclei formation in water is negligible and, in the absence of sites for 

20 heterogeneous nucleation, boiling within a 25 ms pulse is not expected. If nucleation 
sites are present in the liquid remote from the walls, boiling could occur. While it is 
possible that the ZP or local suspended particles may provide sites for nucleation, it is 
unlikely since sharp solid edges are generally more favorable. 

25 Thermal calculations indicate that highly superheated water is briefly formed in 

the focal waist of the IR beam. Heat is conducted mainly radially away from the waist, 
and regions near the beam will be heated by conduction. Superheated water is an 
excellent solvent, and it is not surprising that the region of ZP within the focal waist 
rapidly disappears. 

30 
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After the beam is turned off any vapor bubbles formed would be expected to 
collapse in less than 1 ms, releasing heat to the liquid. If a significant degree of phase 
change had occurred, the temperature at the beeufn core would be held to close to 100°C, 
which would result in lower radial temperatures in the focal plane. The presence of 
superheating could be tested by measuring local temperature excursions. 



Table 1 Peak Central Temperatures 



Power (mW) 


Pulse Duration (ms) 


At End Pulse (°C) 


Steady State (°C) 


100 


1.5 


155.4 


233.9 


100 


3.0 


172.3 


233.9 


100 


4.5 


182.3 


233.9 


100 


6 


189.3 


233.9 


100 


10 / _ 


201.9 


233.9 


ioo 


15.0 


211,7 


233.9 


100 


25 


223.0 


233.9 


80 


15.0 


176.8- 


194.5 


40 


15.0 


106.9 


115.8 



Example 5 : Temperature Check From Probe Beam Steering 

Beam Steering 

In order to confirm the predicted temperature excursion, the change in local 
refractive index around the IR focal waist caused by the temperature excursion was 
estimated. The steering effect caused by the thermal-induced refractive index gradient 
on an orthogonal probe beam was used. Estimates show that a 2° change in ray 
direction would be possible for light passing near the IR heated water column. 

Optical System ' • 

The setup is shown in FIG. 12. The IR beain passes vertically into a 1 x 1 mm 
glass-walled rectangular tube (Vitrocom, Mountain View, NJ) containing distilled H2O. 
The laser is set to pulse at 4 Hz, while beam power and pulse duration are varied. The 
IR beam focus is set at approximately 200 pm above the bottom wall of the rectangular 
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tube. The IR beam focal diameter has already been measured at close to 2a = 6 fim. A 
horizontal probe continuous wave (cw) HeNe laser beam is focused with an f = 120 mm 
lens into the region of the IR beam focus, with the HeNe optic axis orthogonal to the IR 
beam axis. The HeNe probe beani radius to the 1/e^ points at its focal point was 
5 measured by knife intercept as 22.3 ±1.2 |jim. The probe covers the region of interest 
near the IR focus. The HeNe beam is then focused with a cylindrical lens on to a 1 mm 
aperture in front of a Si detector. By scanning the aperture + detector in the direction 
orthogonal to both laser axes, the angle by which the probe beam is refracted (or 
scattered) by the thermal field around the IR beam is determined. The system is used to 
10 measure the maximum angle through which the probe beam is turned at given pulse 

length/power. This is then compared with the angle derived from the predicted thermal 
filed. 

Probe Ray Steering . 

15 The curvature of the probe ray is proportional to (and in the direction oO the 

normalized refractive index gradient V n/n. In the IR beam focal plane the temperature 
field, and consequently the refractive index, in the specimen must be radially synmietric 
around the IR axis. The refractive index will decrease toward the IR axis during the IR 
pulse, and a probe ray traveling in the focal plane of the IR beam will be steered away 

20 from the axis (see FIG. 13). The steering angle \|/ in a centrally symmetric refractive 

20 

field is given by the expression 
Equation 6 

25 

where b is the impact factor, r is radial distance from the IR beam axis, n(r) is the 
refractive index, and R=b/n(R) is the distance of closest approach. 

r 

The refractive index n(T) of superheated water is given by lAPWS^^ from 0 to 
30 150-C. This covers the required range except for the beam center. For higher 




HTS-035DV 



31 

temperatures, the lAPWS data was extrapolated. Using the GE solutions for the thermal 
field to calculate n(r), the expected probe ray steering by the changed refractive index 
during the laser pulse as a function of IR laser power and pulse length was derived. The 
impact factor covers the range 0<b<22 |im due to the large probe beam diameter. 

The maximum beam bending is a measure of the maximum temperature gradient, 
which occurs adjacent to the beam radius and corresponds to the specific temperature 
field. The predicted values for pulse duration 25 and 2.5 ms are shown and maximum 
measured steering angles are given in FIG. 14 for IR power 20-100 mW and pulse 
duration 1-25 ms. The laser power shown in FIG. 14 is that delivered at the focal point, 
allowing for absorption after 200 \im H2O. 

-After the beam is turned off, any vapor bubbles formed would be expected to 
collapse in less than 1 ms, releasing heat to the liquid. Possible phase change is 
neglected in the time scales of interest. Consequently, the predictions derived from 
Equation (1) to compute temperatures in the beam focal plane can be used. 

For higher laser delivered power, the central temperature rises and the thermal 
gradient becomes steeper. The predicted values therefore increase. The measured 
values show a tendency to reach a limiting value at high laser powers in the region 
where the predicted temperature is >2G0°C, implying that the gradient has maximized. 
This is unlikely to reflect phase transition: once phase change nucleates at these super- 
heated temperatures, it would be expected to continue until the medium reaches T = 100. 

If the medium can reach a maximum temperature of only lOO^'C due to phase 

♦ 

change, the predicted refractive index gradient never rises high enough to cause the 
observed beam bending, but flattens out at 0 = 0.6''. Since the observed maximum ray 

bending angle is 9 = 1.8°, the results are consistent with central beam temperature > 
lOO^C. 



I 
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"Melting Point" of Zona Pellucid a 

The diameter of the channels cut in the ZP of three bovine oocytes have been 
measured by pulses of power 100 mW and duration 1.5-25 ms, using a 3 |Lim radius 
beam under standard conditions. The channels cut in the ZP of murine two-cell pre- 
5 embryos have also been measured by the same beam. The channel diameter versus 
pulse duration is shown for both species in FIG. 15. The GE model is applied to 
determine temperature history on the focal plane at various distances from the beam axis 
and the peak temperature calculated at each position is shown in FIG. 15. 

10 The bovine ZP material appears to undergo dissolution (melt) between 140 and 

150^C, the temperature decreasing with longer pulse time. The murine ZP appears to be 
more easily removed and vanishes at about 5-10°C lower temperature than the bovine. 

Descloux and Delacretaz^^ measured the thermal dissolution of mouse ZP at 
15 temperatures 61- 73°C over times up to lOOQ s, and derived a rate constant. This 

measurement is extrapolated to estimate the mouse ZP dissolution temperature for 10 ms 
pulse as 1 1 1°C, lower than the experimental estimate averaged over the ZP from mouse 
pre-embrybs. The minimum temperature that must be attained in the beam center in 
order to ablate human ZP material in the pulse times used is likely to be between 130 
20 andl50°C. 

Example 6: Conditions for Minimizing Collateral Damage 

^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ 

It is assumed that the beam central temperature should be maintained at or above 
about ISO^'C in order to penetrate the ZP within pulse time of 10 ms. At the same time, 

25 the temperature at the surface of the nearest cell should be minimized. The typical 
distance to nearest cell is taken as 20 |bim and use the thermal model above to derive 
temperatures on the focal plane at 20 ^im from the beam axis: in FIG. 16 a graph is 
shown of the peak temperature reached during the pulse as function of pulse duration for 
various beam powers. In the present case, a typical thermal excursion would last on the 

30 order of the pulse duration and have peak aniplitude shown. The cells, therefore, 

undergo rapid heating and cooling of amplitude tens of degrees in a few milliseconds. 
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If thermal damage scales are assumed as a chefnical reaction, then damage would 
be proportional to x • exp(- AQ/RT), where T is pulse duration, T is the niean pulse 
temperature, and AQ is reaction energy. This would imply that the pulse time must be 
5 kept extremely short if the temperature rise is significant. 

Thermal excursion amplitude and duration should be minimized, therefore, pulse 
time is low. Short pulse duration and high beam power give the most efficient system 
for maximizing beam cutting power and minimizing local heating. FIG. 11 indicates 
10 limiting lines corresponding to beam center peak temperatures of 150 and 1 10°C. The 
peak central temperature should exceed the limiting lines for rapid ZP ablation, while 
maintaining minimum peak temperature at 20 |Lim. It is evident that the laser pulse 
duration should be as short as possible, consistent with effective ZP dissolution, i.e., 
preferably below 5 ms. 

15 

■ r 
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